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Abstract 


Having  been  successful  In  our  studies  of  Isolated  (gas-phase)  OH  , NH2  , 
BeH  , CN  , BO  , LI2  , and  NO2  , we  are  now  In  a good  position  to  further  develop 
this  area  by  1)  extending  our  calculations  to  chemically  related  molecules, 
e.g.,  SH  , PH2  , MgH  , etc.,  and  2)  beginning  to  look  Into  the  hydration  (or 
other  solvation)  of  the  species  whose  gas-phase  properties  we  already  know. 

Both  of  these  avenues  of  research  are  feasible.  We  now  have  more  powerful 
computer  programs  (MOLECULE)  which  should  permit  us  to  study  the  heavier  analogs 
mentioned  above.  Moreover,  we  have  also  recently  gained  the  experience  needed 
to  Include  solvation  effects  In  our  studies  of  those  anions  which  we  earlier 
studied  as  Isolated  systems.  The  hydration  may  be  described  either  in  terms  of 
"model"  H2O  molecules  which  are  simply  point  charges  distributed  to  represent 
the  dipole  moments  of  the  surrounding  solvents,  or  In  terms  of  real  ab  Initio 
H2O  molecules.  The  Initio  solvent  studies  will  probably  be  restricted  (by 
cost)  to  cases  in  which  only  one,  two,  or  three  solvent  molecules  surround  the 
anion.  The  "model"  solvent  calculations  might  not  be  sufficient;  we  can  only 
answer  these  questions  by  trying  the  computation,  although  some  existing  work 
on  hydration  by  Morokuma  and  by  ourselves  has  been  encouraging. 

Because  we  have  successfully  carried  out  many  calculations  on  anions 
using  the  method  described  here  In,  a critical  evaluation  of  the  theory  is 
available,  and  therefore,  calculations  on  many  other  chemically  Important  species 
can  now  be  undertaken  with  some  confidence.  In  addition  to  the  above  two 
extensions  of  our  work,  we  are  specifically  interested  in  pushing  on  toward 
obtaining  accurate  electron  detachment  energies,  geometries,  and  charge  densities 
for  the  following  Interesting  ions:  0^  , CO^  , NO^  and  other  anions  which 

play  Important  roles  in  atmospheric  chemistry.  (0^  and  NO2  have  already  been 
characterized).  Now  that  we  have  more  powerful  computer  programs  available. 
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such  studies  are  within  our  range  (with  the  MOLECULE  program  a calculation  on 
NO2  with  30  contracted  Gaussian  orbitals  costs  less  than  $50).  Of  course,  the 
time  and  computer  costs  needed  to  perform  a complete  Investigation  of  all  the 
species  listed  above  may  well  exceed  the  time  Interval  of  this  proposal  and  the 
funds  presently  available  to  us.  However,  the  mere  fact  that  the  geometry 
and  stability  of  0^  remain  as  unsettled  problems  Is  sufficient  to  justify  a 
concentration  of  work  on  this  problem. 

We  also  plan  to  expand  our  research  Into  novel  new  anions  whose  proper- 
ties have  not  been  determined  experimentally.  Our  work  on  the  binding  of  electrons 
to  (closed-shell)  highly  polar  molecules  (LIH,  LIF,  BeO,  LlCl,  NaH,  etc.)  will 
be  extended  to  Include  model  studies  In  which  the  potential  energy  field  of  the 
polar  molecule  Is  replaced  (e.g.  for  a diatomic)  by  that  of  two  point  charges 
(+q)  separated  by  the  experimental  bond  length  (with  Rq  = plus  a "core" 

2 226 

(e.g..  Is  for  LI,  Is  2s  2p  for  Na)  on  the  electropositive  atom.  In  this 
manner  we  hope  to  extend  our  Investigation  to  systems  containing  several  (a 
cluster)  such  polar  molecules.  Our  Initial  results  from  this  kind  of  a model 
calculation  are  very  encouraging. 


I 
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I . Historical  Perspective 

In  1968  Berry^  reviewed  the  experimental  and  theoretical  progress  which 
had  been  made  toward  understanding  the  stabilities  and  bonding  characteristics 
of  small  isolated  (gas-phase)  negative  ions.  In  this  review  Berry  commented 
as  follows: 

’’For  the  theorist,  electron  affinities  and  other  properties  of  negative 
ions  pose  greater  difficulties  than  do  properties  of  neutrals  or  positives, 
insofar  as  electron  correlation  plays  a relatively  larger  part  in  determining 
the  properties  of  a negative  ion  than  it  does  in  other  species.  In  fact, 
electron  affinities  are  frequently  about  the  same  size  as  the  differences 
between  correlation  energies  in  atoms  and  in  the  corresponding  negative  ions." 

As  an  example  of  the  magnitudes  of  electron  correlation  effects,  one  need 

only  consider  the  results  of  our  calculations  on  the  vertical  (R  = 1.718  au) 

electron  detachment  energy  (X  2.  OH  -►  X OH)  of  OH  . Using  an  atomic 

orbital  basis  consisting  of  twenty  Slater-type  orbitals  (STO's),  we  obtained 

a Koopmans'  theorem  approximation  to  the  detachment  energy  equal  to  3,06  eV. 

The  energy  difference  between  two  separate  SCF  calculations  (ASCF),  one  on  OH” 

and  one  on  OH,  carried  out  within  the  same  basis  was  equal  to  - 0.20  eV;  the 

difference  between  Koopmans'  theorem  and  ASCF  represents  the  effects  of  allowing 

the  orbitals  to  relax  upon  removal  of  the  tt  electron.  Finally,  our  best  computed 
2 

energy  difference  , which  contains  effects  of  electron  correlation  through 
third  order,  was  1.76  eV  which  is  in  good  agreement  with  both  Branscomb's 
early  experiments^  and  Lineberger's  more  recent  laser  detachment  results.**  The 
difference  between  1.76  eV  and  the  ASCF  value  of  -0.2  eV  represents  the  effects 
of  electron  correlation.  These  effects  are  indeed  as  large  as  the  entire 
electron  affinity  of  OH;  moreover,  this  result  is  not  atypical. 

Because  the  treatment  of  both  orbital  relaxation  and  electron  correlation 
effects  in  a sufficiently  rigorous  manner  is  an  absolute  necessity  in  any 
reliable  scheme  for  computing  properties  of  anions,  theoretical  progress  toward 
understanding  negative  ions  has  been  made  rather  slowly.  Quite  simply  put. 
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it  is  difficult  to  include  correlation  effects  to  a high  enough  order  to 
guarantee  precision  ± 0.2  eV  in  computed  ion-neutral  energy  differences.  In 
Berry's  review  article,  his  assessment  of  the  state  of  quantum  chemical  research 
on  anions  involved  briefly  mentioning  the  works  of  Pekeris^  on  H”,  Weiss^  on 

Li  , Na  , and  K , dementi  and  co-workers  on  several  atomic  ions,  Sinanoglu 

- - - 9 - 10 

on  C , 0 , and  F , Taylor  and  Hams  on  , Wahl  and  Gilbert  on  halogen 

diatomics,  and  Cade^^  on  OH",  CH",  SiH",  SH~,  and  PH".  Therefore,  in  1968  it 

would  have  been  fair  to  say  that  negative  molecular  ions  were  not  yet  within 

the  class  of  species  which  could  be  conveniently  studied  by  existing  quantum 

chemical  methods.  On  the  other  hand,  the  development  of  modern  laser  technology 

was  making  available  to  the  experimentalist  new  tools  to  use  in  carrying  out  high 

precision  photodetachment  and  photoelectron  spectroscopy  studies  of  gas-phase 

anions.  Thus,  even  in  1968  a great  deal  of  experimental  progress  was  beginning 

to  be  made.  These  experimental  developments  made  a parallel  development  of 

theoretical  methods  and  models  aimed  at  better  understanding  negative  ions  a 

necessary  and  quite  natural  step  in  the  scientific  progress  in  this  area. 

. . 12 

In  1973  Simons  and  Smith  published  an  article  in  which  they  attempted 

13 

to  use  equations  of  motion  (COM)  techniques  to  express  the  vertical  electron 

affinity  (EA)  or  detachment  energy  (DE)  of  a closed-shell  species  in  a manner 

which  treated  orbital  relaxation  and  electron  correlation  through  third  order 

in  perturbation  theory  (the  difference  between  the  coulombic  interaction  and  the 

Hartree-Fock  interaction  being  the  perturbation).  This  developmental  paper 

14-15 

was  followed  by  other  formal  papers  by  Simons,  J«Jrgensen,  Jordan,  and  Chen 
in  which  small  deficiencies  in  the  original  theory  were  corrected  and  connection 
made  with  the  recent  Green's  function  developments  of  Cederbaum  , Pickup  and 
Goscinski^^,  Purvis  and  Ohrn^^,  and  Freed^^.  The  result  of  these  papers  was  a 
method  which  permits  the  direct  calculations  of  EA's  and  DE's,  of  closed-shell 
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species,  which  are  accurate  through  third  order. 

2 20~26 

In  succeeding  publications  ’ ~ the  third  order  EOM  theory  was  applied 

to  studies  of  the  stability  and  bonding  characteristics  of  several  molecular 
anions  (and  cations).  Calculations  of  the  electron  affinities  of  OH,  BeH,  NH^, 

CN,  and  BH  provided  theoretical  support  (to  within  ± 0.2  eV)  for  existing 
experimental  measurements.  Studies  of  the  EA's  of  BO  and  Li^  resulted  in 
theoretical  predictions  for  species  where  good  experimental  data  is  not  available. 
Calculations  of  the  EA's  of  LiF,  LiCJ.,  LiH,  NaH,  and  BeO  have  led  to  predictions 

to  both  the  existence  and  the  stabilities  (with  respect  to  electron  loss)  of 

. . 92 

the  anions  of  these  species.  Jordan  has  also  examined  the  dimer  anion  (LiH)^  ; 
93 

Simons  and  Jordan  have  recently  found  to  be  stable  even  though  Be^  is 

- 27 

unbound.  Very  recently,  the  ion  LiCil  was  observed  by  Lineberger  thereby 

77 

verifying  the  theoretical  prediction  of  Jordan  et  al^  . The  precise  nature  of 

some  of  the  above  calculations  together  with  the  principal  conclusions  of  these 

works  are  treated  later  in  this  proposal. 

In  the  time  since  Berry's  1968  review  article  was  completed  and  for  a few 

years  prior  to  that  date,  several  ab  initio  calculations,  in  addition  to  those 

mentioned  above,  were  performed  on  molecular  anions  which  are  of  chemical 

interest.  These  studies  include  the  following  works:  dementi  (N-  ), 

o 

29  - 30  - - 31  - 32  - 

Lipscomb  (PO  ),  Krauss  (BH^  , 0^  ) Kaufman  (0^  ),  Csmizmadia  (CH^  , NH^')* 

33  - 3U  - 3S 

Pfeiffer  (NO^  ),  Popkie  (C^  ),  McLean'"'’  (0CN~,  SCN"),  Wahl  (012”,  F2'), 

rink^"^  (OH",  NH2",  CHj",  BH^"),  Geller^®  Thulstrup^^  (N0~), 

40  — 41  — 42  ftli 

Peyerimhoff  )♦  Heaton  (NH2  ),  Schaefer  (NO^),  Thulstrup  (C^~), 

Q 5 — — — 

O'Hare  (NF  , NS  , PF  , SF  ) . This  list  is  by  no  means  a complete  tabulation 
of  all  work  done  on  negative  molecular  ions;  it  is  simply  meant  to  indicate 
the  kinds  of  systems  which  have  been  studied  as  well  as  approximately  the  number 
of  calculations  which  have  been  performed  to  date.  Although  it  is  true  that 
a reasonably  large  number  of  ab  initio  calculations  have  been  carried  out 
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for  diatomic  and  small  polyatomic  species,  very  few  of  the  studies  referenced 
above  include  any  electron  correlation  effects.  Most  of  these  calculations 
have  been  done  at  the  SCF-level  within  small-to-moderate  size  bases.  Therefore, 
the  EA's  which  have  been  obtained  in  this  manner  are  probably  not  reliable.  On 
the  other  hand  the  equilibrium  geometries  and  charge  densities  obtained  in  the 
above  SCF-level  works  may  not  be  any  less  accurate  than  the  results  of  analogous 
calculations  on  neutral  species;  electron  correlation  effects  are  not  as  dominant  as 
in  determining  EA's.  Nevertheless,  it  is  our  feeling  that  most  of  the  reliable 
work  on  negative  ions  has  been,  and  will  continue  to  be,  characterized  by  a 
careful  treatment  of  electron  correlation  and  charge  relaxation.  For  this 

n,  the  remainder  of  this  proposal  will  be  restricted  to  discussion  of  results 
proposals  for  studies  which  treat  correlation  in  an  ^ initio  manner, 
ept  in  those  cases  where  simple  chemical  models  are  introduced  to  permit  a 
qualitative  understanding  of  more  complex  systems. 

With  this  brief  survey  of  the  developments  made  since  1968  as  a background, 
let  us  now  turn  to  a more  detailed  discussion  both  of  the  most  recently 
utilized  methods  and  the  results  which  have  been  obtained  with  these  methods. 

In  Sec.  II  we  review  the  foundation  of  the  direct-calculation  approach  of 
refs.  12-19.  Sec.  Ill  contains  a survey  of  results  on  OH  , NH^",  BeH  , CN~, 

BO  , LiF  , LiCA  , LiH  , NaH  , BeO  , NOj  in  which  the  effects  of  orbital 
relaxation  and  correlation  have  been  included.  Specific  attention  is  paid  to 
stabilities  (EA  or  DE),  geometries  (Re,  9^)*  vibrational  frequencies, 
dissociation  energies,  and  charge  densities  of  these  anions.  In  Sec.  IV  we 
review  the  conclusions  which  have  been  reached  thus  far  and  we  make  concrete 
proposals  for  research  in  areas  which  seem  to  show  special  promise  for  future 
development. 
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II . Review  of  the  Theoretical  Methods  Employed 


The  electron  propagator  or  the  one-electron  Green's  function  has  been 

, - 4.-  16,18.20-26,43-46  . ^ ^ ^ , 

used  for  some  time  in  the  study  of  electron  spectroscopy. 

The  advantages  of  using  the  electron  propagator  arise  because  the  transition 

energies  and  the  transition  strengths  are  obtained  directly  as  poles  and 

residues  of  the  propagator,  respectively.  Several  alternative  procedures  for 

decoupling  the  Equation-of-Motion  for  the  electron  propagator  have  been 

developed.  In  this  analysis  we  use  the  superoperator  formalism  of  Goscinski 
47 

and  Lukman  as  the  framework  for  our  development  of  an  electron  propagator 

which  is  consistent  through  third  order.  In  an  alternative  derivation  using 

the  Equation-of-Motion  formalism  of  Rowe  , Simons  and  Smith  attempted  to 

obtain  an  Equation-of-Motion  which  was  consistent  through  third  order.  Purvis 
49 

and  Ohrn  pointed  out  some  deficiencies  in  the  theory  of  Simons  and  Smith; 

these  deficiencies  are  mentioned  again  in  this  section.  We  show  further  how 

the  electron  propagator  can  be  obtained  consistent  through  third  order.  The 

consistency  is  made  more  transparent  by  demonstrating  that  all  second  and  third 

16 

order  self  energy  diagrams  of  Cederbaum  are  included  in  our  formalism. 
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The  definition  of  the  spectral  electron  propagator*'"^  can  be  written 
within  the  superoperator  formalism  as 


G (E)  = (aUEI  - H)“-^la) 


where  I and  H are  the  superoperator  identity  and  hamiltonian  respectively 
and  the  £ are  a set  of  annihilation  operators  a = {a^},  which  are  arranged 

in  a super  row  vector.  The  superoperator  scalar  product  is  defined  in  the 

. . 47  - - 

conventional  fashion  . The  superoperator  resolvent  (El  - H)  can  be 

approximated  via  the  inner  projection  technique  and  the  propagator  then  takes 

the  form 


G (E)  = (ajh)  (hjEI  - H |h)“-'-(h la) 


where  h is  a projection  manifold  which,  if  chosen  to  be  complete  and  ortho- 
normal,  makes  Eqs.  (1)  and  (2)  identical.  The  operator  space 


{hi,h3,h^. . .}  = {a^,ataj^a^,a'^ataj^aj^a^. . . },  k>l,  i>  j , k>l>r... 


spans  the  manifold  h.  We  now  discuss  appropriate  selections  of  h which,  in 
conjunction  with  our  choice  of  the  reference  state,  ensure  that  the  electron 
propagator  is  calculated  correctly  through  third  order  in  the  electronic  inter- 
action. 

. 17  45 

It  IS  well  known  ’ that  the  projection  manifold  h , h„  in  connection 

X ““u 

with  the  Hartree-Fock  (HF)  ground  state  is  able  to  give  the  electron  propagator 
correct  through  second  order  in  the  electronic  interaction.  Our  experience^ 
tells  us  that  second  order  calculation  of  EA's  are  not  sufficiently  precise  to 
be  useful,  (except  in  well-understood  special  cases  which  are  discussed 
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below).  We  demonstrate  how,  using  a correlated  ground  state  and  the  same 
projection  manifold,  we  are  able  to  get  the  electron  propagator  correct  through 
third  order  in  the  electron  interaction. 

The  effect  of  including  h^  in  the  projection  manifold,  where  the  HF  ground 
state  is  used  as  reference  state,  has  been  discussed  by  Tyner  et  and  from 

their  analysis  it  is  clear  that  ^ introduces  terms  which  are  at  least  fourth 
order  in  the  electronic  interaction,  independent  of  the  choice  of  reference 
state.  We  therfore  concentrate  on  using  as  our  projection  manifold  in 

our  search  for  a theory  which  is  consistent  through  third  order. 

As  the  reference  state  in  our  analysis  we  use  a con  olated  wavefunction 
given  by 
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where  the  ^ are  a set  of  HF  creation  operators  and  where  indices  m,n,p,q 
(a,e,6,Y)  refer  to  unoccupied  (occupied)  spin  orbitals  in  the  HF  ground  state 
and  i,j,k,l,r  are  unspecified  spin  orbitals.  We  take  the  correlation 
coefficients  from  Rayleigh-Schrodinger  perturbation  theory 
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where  the  first  number  in  the  bracket  indicates  the  lowest  order  in  the 
electronic  interaction.  The  g^  indicate  HF  orbital  energies  and  the  two 
electron  integral  <mn|a6>  refers  to  the  charge  densities  ma  and  ng  and 

<mn||a6>  = <mn|a6>  - <mn|6a> 


in  our  analysis 
for  convenience 


we  consider  the  projection  manifold  ihj^.h^}  where  the  h 
has  been  redefined  as 


space 


h^  = {a^a^a^  + 


(9) 


The  choice  of  the  subspace  ^ ensures  that  this  space  is  orthogonal  to  h^ 


(10) 


even  for  a correlated  reference  state.  We  also  have  the  following  orthogonality 
relations 


= 4.  ^1^3)  = E (0,2,3...)  (11) 

Using  Eqs.  (10)  and  (11),  Eq,  (2)  can  be  partitioned  into  the  form 

G'^  (E)  = (hjEI  - H|h^)  - (hjH|^)(hg|EI  - H|h^)"^  %|h|11i) 

= Ei  - 4 - B £ 


(12) 
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where  the  matrices  ABC  and  D are  defined  as 


4 = (hjHlUi) 

a = (}^|EI  - Hlh^) 

(13) 


We  will  now  make  an  order  by  order  analysis  of  Eq.  (12)  in  which  we  retain  only 

those  terms  which  are  zeroth,  first,  second,  or  third  order  in  the  electronic 

• • • 12 
interaction.  Since  the  B and  £ matrices  are  at  least  of  first  order  , we  need 

to  consider  only  that  part  of  the  ^ matrix  which  is  zeroth  and  first  order. 

This  constrains  the  indices  in  the  projection  manifold  h_  to  be  of  the  form 

•— O 

a^a  a.  a>6  or  a^a  a m>n  since  operators  such  as  a^a  a lead  to  matrix  elements 
map  amn  mna 

in  the  B matrix  which  are  at  least  of  second  order.  The  resulting  subspace  h. 
is  thus  identical  to  that  used  in  calculating  the  electron  propagator  correct 
through  second  order.  To  calculate  the  electron  propagator  through  third  order 
in  the  electron  repulsions,  we  thus  need  to  obtain  the  £ 3nd  B "'’^trices  through 
second  order,  the  B matrix  through  first  order,  and  the  A matrix  through  third 
order. 

12 

Previous  attempts  to  obtain  the  electron  propagator  correct  through  third 
order  have  used  as  reference  state 


0 = (1  + y k"*"  a^a%  a-)|HF> 

0 aB  m n a 6 ' 

a>B 


A 


(14) 
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theory,.  In  calculating  the  (g.)  and  g,  matrix  elements  correct  through  second 
and  first  order,  respectively,  no  changes  are  obtained  from  considering  the 
higher  correlated  ground  state  in  Eq.  (M).  The  matrix  elements  of  B,  g,  and 
D are  given  by 


B- 

i,am6 


-<im|aB> 


h I <iin|pq>K^P  + I [<1-^1  pa> 

p.q  Y.P^ 


(15) 


B. 

1 ,nam 


<ia  mn>  + 


I 

y6 


<ia 


16y>  k' 


mn 

6y 


I <ip|Yn>K' 

YpL 


mp 

dY 


(16) 


= B (through  second  order) 


(17) 


D 


nam,6pY 


= 0, 


(18) 


~ ^ r, ^ - 0 ) -6  <m6 1 1 pa>  - 6 <nBllqa> 

ndTn^Qop  HQ  dp  nip  m n ci  qn  * pm  ^ 

'^qm  ^ <mnllpq>  + 6^^  <mBijqa>,  (19) 

V.aqB  = -^^6aVYB  ' \b  "'5ql|c»p>  - 6^„<Yq||6p> 

+ ^Ya^^qMBp>  - 6^^<5Y||Ba>  + 6^g<Yq|  |ap>  (20) 

In  the  4 matrix  we  need  to  include  all  terms  up  to  third  order.  The  4 matrix 
elements  obtained  by  using  Eq.  (14)  as  a reference  state  need  to  be  modified 
by  third  order  terms  which  result  from  interaction  between  the  singly  excited 
states  and  the  HE  ground  state.  The  triply  excited  states  which  also  result 
from  a second  order  Rayleigh- Schrodinger  perturbation  calculation  do  not 
introduce  third  order  terms.  We  thus  have  to  add  to  the  4 matrix  elements 
given  in  ref.  [l2]  the  terms  6A.^ 
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y jp  I I i^><6B  mnxmn  pB> 

xa  (e.+L-e  -e  ) 

SB  6 p 6 P m n 


^ I 

pn 


(e,-g  ; 

6 p 

1 (E,+E„-e  -E  , 
0 B m n 

I 

<jp|  |i6><6n|  iBaxaSl  1 

Pn> 

(e.-e  J 

6 p 

1 (e  +Eq-e  -e  ! 

a B p n 

1 

<j«l  Ur 

)XBoi  1 1 6nxpn  | I 

|a6> 

1 (e  +e„-E  -E 

a B p n 

) 

The  A matrix  of  ref.  [l2]  is 


A..  = 6..E.  + y <ik|jl>F,  , 

ID  XD  1 kl 

F = y rP^  - y rP'^  rP^I 

a>6,p  p<q,a 


We  thereby  have  calculated  the  electron  propagator  consistently  through  third 
order. 

A comparison  with  a diagrammatic  perturbation  expansion  of  the  self- 
energy®^  makes  it  further  transparent  that  we  have  really  included  all  terms 
through  third  order  in  our  analysis  of  the  electron  propagator.  In  Fig.  1 we 
have  displayed  the  terms  A3  - A6  as  diagrams  using  the  rules  of  Brandow®^  which 
combine  the  Goldstone  diagrams  with  the  antisymmetrized  vertices  of  Hugenholtz. 
The  diagrams  labelled  A3  - A6  are  identical  to  the  third  order  diagrams  given 
by  Cederbaum  in  which  dots  refer  to  antisymmetrized  vertices.  The  diagrams 

..  49 

A3  - A6  were  shown  by  Purvis  and  Ohrn  to  be  the  only  missing  third  order 


L 


I 


i 
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diagrams  which  evolve  from  a theory  where  the  reference  state  is  given  by 
Eq.  (14).  The  analysis  of  Purvis  and  Ohm  considered  (h,|EI  - Hlh_)  as  two 
matrices;  u which  contains  E and  orbital  energies  (zeroth  order  terms)  and  ^ 
which  contains  the  two  electron  integrals  (first  order  terms).  Making  use  of 

the  identity  (u,  + ^)^  = ^,  + ^)  ^ to  expand  the  inverse,  Purvis  and 

* * 1 

Ohrn  identified  the  diagrams  Cl  - C6  and  D1  - D6  of  Cederbaum  to  originate 

from  the  term  while  the  diagrams  A1  - A2  were  found  to  arise  from  the  4 

matrix  previously  given  in  ref.  [l2].  We  thus  have  accounted  for  all  third 

order  diagrams  which  arise  from  an  expansion  of  the  self  energy.  The  electron 

propagator  calculation  thus  needs  to  use  the  second-order  correlated  ground 

state  given  in  Eq.  (4)  as  a reference  state  and  the  subspace  given  by 

Eq.  (3)  as  tb«  projection  manifold  to  be  correct  through  third  order. 

We  have  shown  how  the  electron  propagator  can  be  obtained  correctly  through 

third  order.  Our  development  stresses  the  fact  that  a complete  treatment  of 

the  inverse  of  the  matrix  is  needed  to  guarantee  that  all  desired  terms  are 

included.  Computational  applications  have  so  far  only  been  carried  out  using 

a diagonal  approximation  to  the  D matrix.  This  situation  is  unsatisfactory 

and  should  be  improved.  A unitary  transformation  which  brings  the  D matrix 

closer  to  diagonal  form  is  related  to  the  theory  of  linear  response  as  discussed 

by  J«5rgensen  and  Purvis  . Using  this  kind  of  procedure  we  would  expect  to  get  I] 

an  approximation  to  the  propagator  which  would  be  nearly  complete  through  j 

third  order.  The  energy  shifted  denominators  which  result  from  evaluating  the  j 

matrix,  correspond  to  the  result  of  summing  certain  classes  of  diagrams  to  j 

i 

infinite  order,  which  implies  that  the  electron  propagator  treatment  has  the'  ; 

i 

computational  advantage  of  expressing  these  summations  in  closed  form.  In  a 
diagrammatic  summation  of  self -energy  diagrams  one  has  to  explicitly  account  ■ 

for  each  energy  shifted  denominator  through  each  order.  We  finally  note  that 

i 

the  side  shifted  diagrams  given  in  Fig.  1 do  not  appear  in  the  original  third  ■ 
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12 

order  theory  of  Simons  and  Smith  for  calculating  ionization  potentials  and 
electron  affinities.  These  diagrams  result  from  using  a more  highly  correlated 
wavefunction  than  the  one  considered  by  Simons  and  Smith  as  the  reference  state. 

With  this  derivation  of  the  needed  third  order  equation  accomplished,  let 
us  now  turn  to  ore  detailed  description  of  selected  results  which  have  been 
obtained  by  maj<;'  ^ use  of  the  propagator  approach.  Recall  that  we  are  limiting 
the  discussion  to  studies  which  have  treated  electron  correlation  effects  in  an 
ab  initio  manner.  The  anions  which  we  have  chosen  to  discuss  in  some  detail  can 
be  divided  into  several  classes.  NH^  , OH  , BeH  , are  hydrides  whose  parents 
have  a half-filled  orbital  to  which  the  "extra"  electron  is  added.  CN~,  BO”, 
and  NOj  are  ions  whose  parents  also  have  half-filled  orbitals.  All  of  the  above 
ions  are  closed-shell  species.  Li^~  and  open-shell  ions  which  are 

formed  by  adding  an  electron  either  to  an  antibonding  (Li^”)  or  bonding  (Be^") 
molecular  orbital  of  the  closed-shell  parent.  LiF~,  LiCl”,  LiH~ , NaH~ , and  BeO~ 
are  each  formed  by  adding  an  electron  to  an  essentially  nonbonding  orbital  of 
the  very  polar  neutral  parent;  in  these  ions  correlation  and  relaxation  are  not 
very  important  because  the  "extra"  electron  resides  in  an  orbital  which  is 
localized  on  the  "back"  end  of  the  electropositive  atom  where  it  encounters 
little  dynamic  interaction  with  the  other  electrons.  It  is  clear  that  our 
studies  of  anions  have  covered  a wide  range  of  chemically  quite  different 
species  in  which  the  factor  (or  factors)  governing  the  binding  of  the  electron 
to  the  parent  include  dipole  forces,  valence  attraction  of  a half-filled 
orbital,  polarization  of  charge,  and  formation  of  new  bonding  orbitals  (Bej"). 


17 


III . Brief  Survey  of  Some  Results 

A • OH  : Our  First  Comparison  With  Experiment 

In  carrying  out  the  calculations  on  OH  described  here  , we  employed  an 

atomic  orbital  basis  consisting  of  Slater-type  functions  whose  orbital  exponents 

were  taken  from  the  bases  of  Cade^^  for  HO" , and  of  Cade  and  Huo^**  for  OH. 

Information  describing  our  basis  and  the  essential  results  of  the  SCF  calcu- 

It*  ^ “ 

lation  on  the  parent  X ^ OH  for  this  basis  are  given  in  table  Al.  Note  that 

the  basis  used  in  this  work  is  not  very  large. 

As  shown  in  table  A2,  the  vertical  detachment  energies  computed  using 

the  basis  of  table  Al  are  within  0.10  eV  of  the  experimental  result  quoted  by 
. 4 

Lineberger  e^  al^  . An  important  observation  which  should  be  made  here  is  that 

the  basis  given  in  table  Al  is  capable  of  yielding  a very  accurate  detachment 

energy.  Our  results  show  that  the  theory  of  molecular  electron  affinities 

and  ionization  potentials  in  ref.  (12)  is  capable  of  yielding  the  vertical 

Ir-t  - 

electi’on  detachment  energy  of  X } OH  to  within  0.1  eV.  It  has  also  been 

demonstrated  that  a highly  accurate  description  of  the  core  orbital  of  OH"  is 

not  essential  to  an  accurate  calculation  of  the  valence  electron  detachment 

energy.  Finally,  an  investigation  of  the  roles  of  orbital  relaxation  and 

correlation  energy  change  in  determining  the  ion-molecule  energy  difference 

2 

has  led  to  the  conclusion  that  both  of  these  effects  must  be  treated  properly 
in  any  study  of  negative  molecular  ions  unless  one  knows  that  the  "extra" 
electron  is  essentially  uncorrelated  (perhaps  by  spatial  localization)  with  the 
other  electrons.  As  discussed  in  part  E of  this  section,  such  is  the  case  for 
the  family  of  anions  formed  by  adding  an  electron  to  a highly  polar  closed-shell 


.-.x. 


molecule. 


(2.4385)  -.0014 
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Table  A2 


Surmary  of  Detachment  Energies  for  OH" 


) 


Detachment  Energy  (eV) 

i 

Source 

1.773 

1 

EOM^ 

1 

! 1.825  ± 0.002 

Hotop,  Patterson,  Lineberger^ 

1.83 

1 ... 

Bransccffnb^ 

a.  Hotop,  Patterson,  Lineberger,  J.  Chem.  Phys., 

b.  L.  M.  Branscomb,  Phys.  Rev.,  1^,  11  (1966). 

M,  1806  (1974). 

i 

J 

'I 
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B . BeH  : An  Experimentally  Understudied  Anion 

An  initial  basis  set  for  the  closed-shell  BeH"  consisting  of  20 

Slater-type  orbitals  (STO's)  was  adapted  from  the  optimized  BeH  basis  set 

54 

reported  by  Cade  and  Huo.  To  accomodate  the  extra  electron  correlation,  2p 
functions  and  diffuse  s and  2p^  functions  were  added  to  the  "sigma  only" 

BeH  basis  set  to  replace  functions  contributing  nominally  to  the  description  of 
the  occupied  BeH  molecular  orbitals. 

The  orbital  exponents  of  the  four  BeH"  STO's  in  the  original  basis  set 
having  the  largest  expansion  coefficients  in  the  Bo  highest  occupied  molecular 
orbital  (HOMO)  were  optimized  at  the  initially  calculated  BeH~  equilibrium 
internuclear  distance  of  2.660  au.  The  initial  and  optimized  BeH"  basis  sets 
and  expansion  coefficients  for  occupied  molecular  orbitals  are  listed  in 
Table  Bl.  Basis  functions  which  were  also  used  in  the  Cade  and  Huo  BeH 
basis  set  have  been  marked  with  an  asterisk.  From  this  table  we  observe  that 
optimization  of  the  BeH  basis  set  cuased  dramatic  increases  in  the  importance 
of  the  diffuse  2sBe  and  IsH  basis  functions  describing  the  3o  HOMO.  SCF 
calculations  for  Egj,j^-(R)  were  executed  on  the  University  of  Utah  Univac  1108 
computer  using  a modified  version  of  the  Harris  DIATOM  program.  Execution  time 
for  each  run  was  approximately  4.50  minutes.  Vertical  ionization  energies  of 
BeH  , were  calculated  by  our  third  order  equations-of-motion  program. 

Execution  time  for  ea_h  ionization  energy  calculation  on  the  Univac  1108  was 
approximately  25  seconds.  BeH  energies  were  calculated  by  adding  the  vertical 
ionization  energy  of  BeH  to  the  BeH"  energy, 

* ‘'’BeH-"" 


(23) 


1.4500 
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Approximate  vibrational  force  constants  (Table  B2)  and  fundamental 
vibrational  frequencies  for  the  BeH  and  BeH"  systems  were  obtained  by  fitting 
a least  squares  quadratic  polynomial  to  each  potential  energy  curve.  Agreement 
of  these  calculations  with  the  experimental  BeH  values  was  quite  good.  For 
the  original  basis  set,  (calc.)  _ (exp.)  +44.6  cm"^,  a +2.17% 

deviation,  and  for  the  optimized  basis  set  the  deviation  was  -60.0  cm”^  or 

BcH 

-2.91%.  For  comparison,  the  value  for  reported  by  Cade  and  Huo  deviated 

from  the  experimental  value  by  +95.0  cm“^  or  +4.66%.  It  is  reasonable  to  assume 
our  calculated  fundamental  vibrational  frequencies  for  BeH~  are  at  least  this 


accurate . 

In  reporting  Feldmann's  value  for  IP_  „-(Thermo.)  of  0.74  eV  in  Table  B2, 
we  have  subtracted  0.01  eV  (zero-point  corrections)  from  the  BeH  photodetachment 
enei'gy  value  of  0.75  eV.  Calculated  thermodynamic  ionization  potentials  for 
BeH  deviate  from  the  experimental  value  by  [-0.01  eV, \ [+0.03  eV, 
+4.1%]^°^’ \ Vertical  ionization  energies  and  electron  affinities  differ  no 
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C . CN  and  BO  ; A Well-Known  Ligand  and  its  Isoelectronic  Partner 

The  starting  point  for  the  construction  of  the  double  zeta  basis  sets 
used  in  these  calculations  was  Roetti  and  dementi's  excellent  set  of 
double  zeta  functions  for  the  component  atoms.  To  better  describe  the  charge 
distribution  in  the  resultant  negative  molecular  ions,  the  orbital  exponents 
of  the  functions  with  large  expansion  coefficients  in  the  higliest  occupied 
molecular  orbital  (HOMO)  of  CN  and  of  BO  were  varied  to  maximize  the 
ionization  energies.  The  greatest  changes  in  each  case  were  produced  by 
modifying  the  exponents  of  the  2s  and  2pa  functions  on  the  less  electronegative 
atom  in  the  ion. 

Since  little  is  known  about  the  geometry  of  B0~,  its  basis  set  was 
optimized  at  R = 2.278  bohr,  which  corresponds  to  the  equilibrium  separation 
of  the  neutral  BO  molecule.  The  total  energy  before  optimization  was  -99.550 
hartree  and  the  corresponding  \ \ ionization  potential  was  2.16  eV. 

The  final  basis  gave  a slightly  lower  SCF  energy,  -99.554  hartree,  and  a 
vertical  ionization  potential  of  2.81  eV. 

The  CN  basis  set  optimizations  were  performed  , at  R = 2.2  bohr,  the 
internuclear  separation  corresponding  to  'the  lowest  energy  found  after  a 
few  preliminary  SCF  calculations  on  CN  with  the  starting  basis.  The  initial 
energy  at  2.2  bohr  was  -92.2634  hartree,  and  the  vertical  ionization  energy 
for  the  starting  basis  was  3.04  eV.  After  optimization  of  2s  and  2p(j  functions 
on  carbon  and  2s  functions  on  nitrogen,  the  total  energy  was  -92.2645  hartree 
and  the  ) -*■  X ] vertical  ionization  energy  was  3.69  eV.  The  diffuse 

2pTi  function  on  nitrogen  was  varied  slightly  in  an  attempt  to  obtain  a 

Ir  2 . . . . 2 

reasonable  I ionization  potential.  However,  the  it  state  of  CN  is  not 

expected  to  be  well  described  in  our  basis  since  the  optimization  of  the  pi 
functions  was  npt  extensive,  The  final  optimized  CN“  and  B0“  basis  sets 
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shown  in  Table  Cl  were  used  to  compute  ion-molecule  energy  differences  at  • 
several  internuclear  separations;  these  differences  were  added  to  the  total 
SCF  energies  at  corresponding  R values  of  the  negative  ions  to  generate 
SCF-level  potential  curves  for  BO  and  CN. 

Examination  of  the  ground  state  potential  curves  for  BO,  B0~,  CN  and 
CN”  allows  us  to  determine  the  adiabatic  electron  affinities  of  BO  and  CN, 
which  we  can  compare  with  existing  experimental  data.  Our  predicted  electron 
affinity  for  BO,  2.79  ± .2  eV,  is  within  the  range  of  experimental  estimates 
which  vary  from  2.4  eV  to  3.1  ± .1  eV.  Experimental  studies  of  CN  have 
yielded  more  precise  results.  Chupka,  et  al^^,  have  reported  an  electron 
affinity  of  3.82  ± .02  eV,  which  was  obtained  from  photodissociation  measure- 
ments on  HCN.  Our  calculated  electron  affinity,  3.70  ± .2  eV,  is  in  good 
agreement  with  this  value. 

The  results  of  our  EOM  calculations  and  Koopmans ' theorem  estimates 

for  vertical  ionization  potentials  of  CN"  and  BO”  are  shown  in  Table  C2  for 

selected  internuclear  separations.  The  Koopmans'  theorem  values  deviate 

1 + 2 '*■ 

considerably  from  EOM  results  for  the  I I ionization  of  CN  , less  so 
1 + 2 

for  the  E E BO  ionization,  and  are  actually  very  close  to  EOM  results 
Iv^  2 . . . 

for  the  i ^ ionization  of  CN  . 

The  large  difference  between  Koopmans'  theorem  and  the  equations-of- 
motion  result  for  the  ^ ^ ionization  of  CN  indicates  that  orbital 

relaxation  and  electron  correlation  effects  are  important  in  the  process  of 
removing  an  electron  from  the  5a  orbital  of  CN  . Analogously,  relaxation  and 
correlation  effects  appear  to  be  less  important  in  the  ionization  of  BO". 


THE  FINDINGS  IN  THIS  REPORT  ARE  NOT  TO  BE  CONSTRUED  AS  AN  OFFICIAL  DEPARTMENT 
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